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We demonstrate that charge-qubit cluster state generation by capacitive coupling is anisotropic. 
Specifically, horizontal vs vertical nearest-neighbor inter-qubit coupling differs in a rectangular lat- 
tice. We show how to ameliorate this anisotropy by applying potential biases to the array of double 
dots. 
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I. INTRODUCTION 

One-way quantum computing is a particularly attrac- 
tive model for quantum circuits because global entangle- 
ment is accomplished in a single step, and then all sub- 
sequent quantum computation is effected simply by se- 
quential feedback-controlled single-qubit measurements^. 
The globally entangled state that serves as a universal 
substrate for quantum computation is known as a clus- 
ter state and was originally proposed for optical lattices^ 
and subsequently demonstrated with photons^. Efficient 
quantum circuits for one-way quantum computing was 
proposed for solid state devices^. 

Solid-state charge-qubit cluster states offer the excit- 
ing prospect of one-way quantum computing with semi- 
conductors^. Here we show that proposals for peri- 
odic generation of charge-qubit cluster states involving 
double-dot (henceforth 'ddot' as this term emphasizes the 
single-entity nature of the ddot structure) charge qubits 
are complicated by an overlooked inter-qubit coupling 
asymmetry in two dimensions. We remedy this compli- 
cation by showing that the original proposals^ can be 
recovered simply by applying a potential field bias. 

We proceed first by establishing the second-quantized 
Hamiltonian description for the array of quantum dots 
and then showing how the Hamiltonian can be simplified 
to a first-quantized Hamiltonian over ddot charge qubits. 
In the slow tunncling-ratc regime, we show that the first- 
quantized Hamiltonian is well approximated by the ubiq- 
uitous Ising-like Hamiltonian. This Hamiltonian consid- 
erably simplifies the dynamical description and shows 
that the necessary bias of the ddot charge qubits to gen- 
erate cluster states is determined by the number of ddot 
neighbors. Thus a global bias of a large structure will 
lead to periodic evolution of excellent approximations to 
cluster states. Furthermore, by applying different biases 
only to the ddots on the boundary (and not to the ddots 
within), the Hamiltonian induces evolution to the ideal 
cluster state. 

Our aim here is to remedy the deficiency of anisotropies 
in the evolution of ddot charge-qubits. We show that this 
problem can be nearly remedied by a global field bias and 
completely remedied by a global field bias with corrective 
biases applied to the ddots at the boundary. This method 




FIG. 1: (Color online.) (a) A two-dimensional cluster state. 
Each grey ball represents a charge qubit, and lines connect 
nearest neighbors in the horizontal and vertical directions. 
These lines correspond to coupling by controUed-Z operations, 
(b) The qubit represented by a grey ball is expanded to a 
double-dot structure. The logical qubit state |0) is depicted as 
the double-dot structure with and extra charge in the left (L) 
dot. (c) The double-dot structure is modeled as a double- 
well potential. The blue (lower) pair of lines corresponds to 
the energy of the symmetric state for the sharing of the excess 
charge between the two wells, and the red (upper) pair of lines 
corresponds to the energy of the antisymmetric shared-charge 
state. 



for correcting anisotropy is examined numerically for the 
case of 40nm GaAs ddots in a two-dimensional lattice 
formation with a being the distance between two sites of 
the ddot, and dy the distance between both left and 
right sites of two nearest-neighbor ddots in x direction 
and y direction in Fig.[lja), respectively. 



II. CHARGE-QUBIT CLUSTER STATE 

The charge-qubit cluster state (CQCS) in two dimen- 
sions is depicted in Fig. [TJ The standard depiction of this 
cluster state is shown in Fig. [Ija) as a periodic rectan- 
gular lattice of qubits connected to nearest neighbors by 
solid lines. Each qubit state is in H2 = span{|0) , |1)}, 
with |0) the logical zero state and |1) the logical one state. 
Cluster state generation proceeds first by globally trans- 
forming every qubit from the state |0) to the state |+) 
where |±) (|0) ± |1)) /V2. We refer to {|0) , |1)} as 
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the 'standard basis' and {|+) , |— )} as the 'dual basis'. 

After all N qubits in the cluster state are prepared in 
the 1+)®^ state, nearest- neighbor qubits then interact 
via the two-qubit controlled- Z operations, denoted CZ. 
Here Z is the Pauli 'phase' operator. The other Pauli 
operators are the 'flip' operator X, the 'flip+phasc' op- 
erator Y = XZ, and the identity operator 1. This op- 
eration is represented in the two-qubit standard basis 
{|00),|01),|10),|11)} as CZ = diag{f,f,f,-f}. The 
unitary operation CZ is independent of whether the line 
in Fig. [TJa) is horizontal or vertical. 

Fig. [Ijb) shows that the horizontal vs vertical sym- 
metry is in fact broken by the coupling axis of the 
charge qubit, which is represented in Fig. [Ijb) as an ex- 
cess charge in the left or right quantum dot. Although 
the charge-qubit coupling axis could be aligned indepen- 
dently of the orientation of the overall qubit lattice, we 
will treat the case that the charge-qubit coupling axis is 
in the x-direction. This is a physically reasonable case, 
and extending to the case of arbitrary alignment is in- 
volved but not difficult. 

The charge qubit can be created as a semiconductor 
ddot structure^. Other alternatives exist such as the su- 
perconducting charge qubit^ or a pair of dangling bonds 
on a surface^. In any case, the logical states typically cor- 
respond either to the left- and right-well occupancy by 
the excess charge or, alternatively, to the cases of sym- 
metric or antisymmetric charge states between the two 
dots of a ddot charge qubit. 

For coherently evolving charge qubits, Schrodinger's 
equation can be used to describe the dynamics, and 
the potential in Schrodinger's equation is depicted in 
Fig. [TJc). Here we treat the standard basis as corre- 
sponding to left- and right-occupancy; the dual basis then 
corresponds to the symmetric and antisymmetric charge- 
occupancy states. 

The quantum dots are engineered so that each poten- 
tial well has only one bound energy state for the excess 
electron. Due to the Pauli exclusion principle, the num- 
ber of electrons in each well is either zero or one or else 
two electrons with opposite spins. The case of two excess 
electrons in one double-dot structure should be energet- 
ically forbidden by Coulomb repulsion between the two 
electrons to preserve the integrity of the charge qubit. 



III. MODELING THE DYNAMICS 

The goal is to have one excess electron per closely- 
spaced quantum dot pair, but the general picture is that 
each quantum dot can have one or two electrons. The 
restriction of one excess electron must emerge as an ener- 
getically favorable configuration rather than be imposed 
by fiat. The full second-quantized description of the elec- 
trons in the array of quantum dots is given by the ex- 
tended Hubbard model (EHM). 



A. The extended Hubbard model 

The EHM applies to an array of quantum dots whose 
locations in a two-dimensional array are denoted by lat- 
tice coordinates. For Cij the annihilation operator at 
dot site clj the conjugate creation operator, and 

fiij = cljCij the number operator, the dynamics of the 
charge-qubit cluster state can conveniently be described 
by the extended Hubbard HamiltonianiS (As spin is con- 
served, we can, without loss of generality, assume fixed 
spin and ignore this degree of freedom) 

H=Y,En^j+V (1) 
+ 51 Wij^i:yh,jh,,j, ~T,j^,,j, {c\^c,>j: +\v.c^ . 

Here E is the effective on-site energy for each site (i, j), 
which can vary due to local field effects. Tij^i>ji 
is the coherent tunneling rate between sites 
and Wij^'j' is the Coulomb repulsion energy be- 

tween sites and Finally, for Anij^i'j' :~ fiij — 

fiiiji the number-difference operator between sites 
and the potential bias operator is 

V =^ Kj,i'i'Anij-i-jv (2) 

with Vij,i'j' the inter-site ^ [i' ij') potential differ- 
ence. 

In fact Eq. ([T]) describes not just nearest-neighbor in- 
teractions but interactions between all dots with all other 
dots, where the inter-dot couplings Tij^iiji and Wij^i'j' 
arc suitably chosen. For charge qubits corresponding to 
closely spaced dot pairs, Tij,i/ji can be neglected for all 
but the ddot of a given charge qubit. Also Wij,i'j' is only 
significant between charge qubits. The intcrdot (possibly 
screened) Coulomb repulsion is neglected for the ddots 
of a charge qubit because Wtj^ij is sufficiently large to 
prevent both dots from being simultaneously excessively 
charged. 

B. Single-qubit gates 

Eq. ([1]) is a second-quantized Hamiltonian. To bridge 
this Hamiltonian over to the multi-qubit description, we 
restrict the Hilbcrt space, upon which the first-quantized 
version of the Hamiltonian acts, to the case of a single 
excess electron in each double well. Note here that a 
dot has lattice coordinates expressed here as («i,ji) and 
dot 2 is at {12,32)- 

For an array with close proximity between dots of the 
ddot pair, the resultant charge-qubit ddot pair can be 
treated as a point-like object in the quantum computing 
architecture. This ddot charge qubit has a point-like co- 
ordinate designated by (m, n) where the change of font is 
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used to indicate point-like ddot coordinates rather than 
the coordinates of a particular quantum dot. Thus, we 
use to designate the location of a quantum dot in a 
two-dimensional array and (m, n) to denote the location 
of a point-like ddot charge qubit. 

Assuming two dots of each ddot pair share one elec- 
tron, the state od ddot pair is a superposition of two 
basis states: {\L) , |-R)}. Here L (R) indicates that the 
electron is in the left (right) dot. In this basis, H in 
Eq. (m for one ddot projected onto Hmn is 

^mn = Ei\L) {L\ + \R) + V{\L) {L\ - \R) {R\ 

+ Ti\L){R\+h.c.))=(^^ + ^ (3) 

for V ^ Vl — Vr (the relative energy between the left 
and right dots of a ddot pair), and T is the flip rate 
corresponding to the tunneling rate between the two dots 
of a ddot charge qubit. 

This Hamiltonian can be conveniently rewritten as 
a linear combination of three types of quantum gates. 
These gates are the identity 1, X = |0) (1| -(- 11) (0|, phase 
gate Z = |0) (0| - |1) (1| with |0) \L) and |1) := \R). 
With these simplifications, the Hamiltonian ([T]) can be 
projected into the qubit space. This projection becomes 
clear by studying the single-qubit case comprising one 
ddot pair. 

In the standard basis the ddot single-qubit Hamilto- 
nian is Hmn = E1+ TX + VZ at site (m, n). Here E is 
an energy term for the qubit. The bias V can be con- 
trolled by applying an electric field potential across the 
ddot pair. 

In order to connect our mathematical expressions to 
a typical experimental setting, we consider a GaAs ddot 
with a single-dot diameter of 40nroi^. A typical experi- 
mental parameter range for tunneling is T ss — lO^eV: 
here we choose T = 0.1/xeV (160MHz). By tuning the 
electric field to ^ = 0, the evolution of the Hamiltonian 
Hmn effectively implements bit flips via the X operator, 
and the resultant tunneling or flipping rate is 160MHz. 
By tuning the electric potential bias to ^ T, e.g. 
V = lOyueV, the dynamics is dominated by phase flip- 
ping at a rate of 16 GHz. 

C. Two-qubit gates 

Now let us consider the two-qubit gate such as the CZ 
gate, which can be implemented via Coulomb interac- 
tion between two nearest-neighbor charge qubits shown 
below. The Hamiltonian for two nearest-neighbor charge 
qubits located at (m, n) and at (m',n') with (possibly 
screened) Coulomb interaction is 

-ffmnm'n' = '^Et + TXmn + KniT^mn + TXm'n' (4) 

1 

+ Vm 

mnm'n' 

{ik\. 

l,k=0 



The last term describes the Coulomb interaction en- 
ergy in the qubit basis between the two nearest-neighbor 
charge qubits located at (m, n) and at (m',n') to with 
sums over I and k representing the two charge qubits are 
in the state \l) and |fc). 

The coefficients Qk SuTC the inter-site Coulomb interac- 
tion strengths between the same or opposite sites of the 
two charge qubits (m, n) and (m', n'). For the case I = fc, 
i.e. the two charge qubits are in the same states and both 
electrons in the left or right dots, we have c^oo = ^ll^ " 
and ?ii = Wfl^'"'"'. For the other case I ^ k, i.e. the 
two charge qubits are in different states and the two 
electrons in different dots, we have c^oi = ^lb^ " a-^^d 

In the 2D lattice shown in Fig. [TJa), for two nearest- 
neighbor charge qubits there are two cases: both qubits 
located in the same column (along the x direction) and 
in the same row (along the y direction). We use the su- 
perscripts X and y to distinguish these two cases. For the 
two nearest-neighbor charge qubits in the same column 
in which the structure is symmetric, c^qq = c^f-^ = Vq/dx, 

and (Toi = ^10 = Vq/V^I + for Vq = eV(47re) with 
e the applicable dielectric constant. In contrast, for the 
two nearest-neighbor charge qubits in the same row in 
which the structure is asymmetric, c^oo = = ^q/rfyj 
'^01 = y^/idy + a), and ^fg = VQ/(dy - a). 

For the GaAs ddot system considered in the previous 
subsection, we have Vq = 1.75 x lO^^^N-m^. With the 
experimental parameters a = 400nm, dx = 5.5/xm and 
dy = 5.85^roi^, we numerically estimate the coefficients 
of the interaction terms to obtain <^gQ = = 20.0^eV, 
^01 = ^fo = 19.8^eV, = 4^ = 18.7//eV, ^q^^ = 17.5//eV 
and 4q = 20.1^eV. 

IV. GENERATING CLUSTER STATES 

In the one-way quantum computing model, the two- 
dimensional cluster state is a highly entangled multi- 
qubit state and processed by performing sequences of 
adaptive single-qubit measurements, thereby realizing ar- 
bitrary quantum computations. The two-dimensional 
cluster state serves as a universal resource for one-way 
quantum computing, in the sense that any multi-qubit 
state can be prepared by performing sequences of local 
operations on a sufficiently large two-dimensional cluster 
stateii. 

In previous proposals^, charge qubits are treated as 
being symmetrically coupled, which is an appropriate 
strategy for the one-dimensional case but not at all for 
the two-dimensional case. Here we show that, by apply- 
ing local external electric fields, we can generate two- 
dimensional cluster states without the requirement of 
symmetry of charge qubits. 

In the two-dimensional case shown in Fig. [Ija), we 
obtain a more general Hamiltonian case with external 
electric fields Vmn applied on each qubit (m, n) located in 
the m*'^ row and the n"' column of the ddot charge qubit 
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array: 



H2D =N'E 



N 



+ 2*^+-"-™" lm,n+l + ^''-^mn ® ^m,n+l 
1 



-(A4 + A^^)Z^„ 
i(A4-A^^)Z^,„+i 



(5) 



For the GaAs ddot system considered in the previ- 
ous section, the energy offsets are = 39.7/ieV, = 
37.5/ieV, A? = -2.6/^eV, and <r = lO.O^eV. Note that A<; 
is comparable to the bias Vmn hence cannot be neglected. 

The Hamiltonian for the two-dimensional array of 
charge qubits can be simplified as (neglecting identical 
terms such as N'^E and other such terms) 

N 

^2D X! [^^mn + KniT^mn (7) 
m,n— 1 

+ "^(^mn ® Zxn+l.n + ^mn ® ^m,n+l)] , 

where 



where in order to rewrite the interaction terms in the 



Pauli operator Z, we introduce 2^^^^^ = <;qq^^ 



''11 ^ 



4^'^ + = ,1^^ - and A<r!(«) = ^i^^ - 

'■10 ■ 

For the two nearest-neighbor charge qubits in the same 
column in which the structure is symmetric, we have ~ 

VQld^ + Fq/VdlT^, e = l^q/rf. and A<j| = 0. In 
contrast, for the two nearest-neighbor charge qubits in 
the same row in which the structure is asymmetric, we 
have 
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dy 



2{dy + a) 2{dy - a) ' 



(6) 



A4 = 0, and A<;^ = A<r = -2VQa/{dl - a^). By choos- 
ing the proper distances dx and dy between the nearest- 
neighbor charge qubits, we have <j^/2 = <j^/2 = <j. 



f Vmn, m= l,...,iV;n = 2,...,7V-l 

^mn = < Kni - A^/2, m = 1 , . . . , iV; n = 1 
[ VmN + Ai/2, m^l,...,N;n^N 

(8) 

are the modified energy offsets for each ddot pair at 
site (m, n). 

For the ddot pairs on the two edge columns, the energy 
offsets are modified due to the asymmetry of the struc- 
ture. In contrast, for those in the middle columns, the 
energy offsets remain and arc caused by the applied elec- 
tric fields because the extra energy offsets are canceled 
out due to the structure. 

V. APPROXIMATING ISING-LIKE DYNAMICS 

We apply a canonical transformation for a global basis 
change on the Hamiltonian shown in Eq. 



Hcs = exp 



I m,n=l J I m,n=l-^'^mn J 



= Hi, + H, 



und ; 



(9) 



where 



of the cluster state) and 



AT 



-^^Is — [£'mn^mn + ^(^mn<8'^m+l,n + ^mn<8)^m,n+l)] 

m.n— 1 

(10) 

for Emn = Knn + ^^^/^n an Isiiig-likc Hamiltonian (we 
need not only the interaction term such as Zmn'^Zm+i,n+ 
Zmn ® Zm,n+i but also the term Z^n for the generation 



Hund 

N 



v 

mn=l "1" 



(11) 



^ -tr, Zm+i,n '8) A^mn ^ TT, Z^n $5 Xm+i,n 



j,2 rj. 

Amn ® Am+i,n + TJ^ — ■Z^m,n+i ® X^in 
+i,n ^vnn 

•2 



^ mn ^ m+i,n 



T 

T7} Zmn ® Am,n+i ^ -TTi — Amn ® Ant,n+i]. 



The Hamiltonian -ffund is an undesirable interaction for 
the generation of cluster states. In the slow-tunneling 
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regime 

T«k|<|Knn| (12) 

combined with T < |V;^„| derived from Eqs. © and ([T^ . 
we obtain the approximate Ising-hke Hamihonian His. If 
T <C \V^ „\ is satisfied, the coefficients of ffund are small 
enough so that the unwanted interaction term i?und can 
be neglected and Emn ~ Krm- Thus, the control term V 
is now embedded within the term i?, which incorporates 
the modified control bias V' and the tunneling rate T . 

A. Periodically generating a cluster state 



TABLE I: For a ddot system in GaAs, with the choice of 
fc = fc' = for simplicity and experimental parameters Vq = 
1.75 X 10"^^N-m^, a = 400nm, = 5.5^m and dy = 5.85^tm, 
we can estimate the potential energy offset for each ddot as 
follows. 



m 


n 


K.n(MeV) 


l,iV 


1 


A<,/2 - 2? = -21.3 


1,N 


N 


-Ac/2 - 2? = -18.7 


l,N 


2,...,7V-1 


-3? = -30.0 


2,...,N-l 


1 


Ac/2 -3? = -31.3 


2,...,N-l 


1 


-Ac/2 - 3c = -28.7 


2,...,N~l 


2,...,iV-l 


-4c = -40.0 



We can periodically generate a large particle cluster 
state 



N 



(13) 



in this tilted (i.e. biased) frame by applying the unitary 
operation exp|— ii?ist| on the initial state l^');^; after 
a time i, if and only if both ci = f + 2fc7r and Emnt ~ 
^§(^^1111^ + ^mn) + 2fc'7r are satisfied for k and k' arbitrary 
integers. Here Vma^ is the number of qubits connected to 
the qubit (m, n) in the x (y) axis shown in Fig. [TJa). 
The two constraints lead to the relation 



Err 



1 



(14) 



Consider a two-dimensional structure shown in Fig. HJa) 
we have 



Ell 
and 

Em 



EiN = Eni = Enn = (8A:' - 2)c/(8fc + 1) (15) 



8fe'-4 



8fe+l ^' 
8fc'-3 
8fe+l ^' 



m, n = 2, . . . 

m = l,iV;n 



= 2, 



..,7V-1, (16) 



8fc'-3 
8fc+l 



c, m 



2,...,A^- l;n= l,iV. 



Equation ([T5|) represents the effective energy offsets of 
the ddot pairs at the four corners, which have two connec- 
tions along the x and y axes respectively; i.e. Vx + t^y = 2. 
The similar principle applies to the first expression of 
Eq. ((T6| . which shows the effective energy offsets of the 
charge qubits that interact with a total of four nearest- 
neighbor qubits: i.e. Vx -\- = 4. The second and third 
expressions of Eq. show that the effective potential 
differences of the ddot pairs on the boundaries with three 
connections; i.e. Vx + Vy = 3. 

From the above reasoning we see that a cluster state 
can be generated for a potential energy offset V^n of the 
(m, n)'^ ddot qubit as 



i^mn, m = l,...,iV;n = 2,.. 

Knn « <( -Bn^i + Ac/2, m = 1, . . . , iV; n = 1, 
^;n^JV-Ac/2, m = l,...,iV;n = A^. 



(17) 




(b) 

















... "'"^ 






FIG. 2: (Color online.) (a) The electric fields with different 
strengths are applied on the ddot pairs in a periodic structure 
isolated by the dashed and solid lines with different colors and 
shapes: (8fc' - 2)c/(8fc + I) + Ac/2 (yellow square); (8A:' - 

2) c/(8fc + 1) - Ac/2 (yellow circle); (8fc' - 4)c/(8/c + 1) (black 
dashed box); (8fc' - 3)c/(8fc -I- 1) (red dashed circle); (8fc' - 

3) c/(8A; + l) + Ac/2 (blue dashed circle); (8A;' - 3)c/(8fc + 1) - 
Ac/2 (green dashed circle), (b) By choosing k = k' = Q for 
simplicity, we show an example of how to apply electric fields 
on the ddot pairs for generating a two-dimensional cluster 
state, (c) The periodic structure of a two-dimensional array 
of ddot pairs for an alternative method to generate a two- 
dimensional cluster state by simply applying a global electric 
field —3c on the first and last rows, and —4c for the remaining 
rows, respectively. 



Hence generating two-dimensional cluster states can be 
achieved by applying local electric fields to set the energy 
offsets as above. In our structure shown in Figs.[2ja, b), 
there are six choices of electric field biases in total with 
different strengths applied to each ddot pair described in 
Table 1. 
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B. Cylindrical cluster state 

By changing the structure of the array of charge qubits 
to that shown in Fig. [5Jc), we can generate a large two- 
dimensional cluster state by simply applying global elec- 
tric fields instead of local electric fields, resulting in an 
important simplification to the technical challenges of 
generating cluster states. 

From our analysis above we observe that the choice of 
external electric field strength depends on two factors: 
the number of connections to nearest neighbors + i^y 
and the extra energy offset due to asymmetry of the two- 
dimensional structure. The periodic boundary condition 
n -I- TV = n dismisses the second factor, and the extra en- 
ergy offsets of the ddot pairs located on two edge columns 
cancel each other. 

The periodic boundary condition also partially dimin- 
ishes the first factor and makes the connections of all ddot 
pairs except for those located on the first and last rows 
to be + i^y = 4. For those ddots located on the first 
and last rows, Vx + i^y = 3, and the global electric fields 
on the two rows are set to be {8k' — 3)<r/(8fc+ 1) while all 
the others are set to be (Sfc' — 4)<r/(8fc 4-1). For simplic- 
ity, choosing k ~ k' = 0, to generate a two-dimensional 
cluster state, one only needs to apply an external global 
electric field — 3^ on the first and last rows, and — 4<r for 
the remaining rows, respectively. 



C. Validity of Ising-like evolution to a cluster state 



Our approach is valid if the contribution of the term 
fl"und in Eq- (HH) to the evolution is neghgible. The fi- 
delity of the pure two-dimensional cluster state is 



T = 1/icV and c = lO^ieV, A^^ 
F = 0.9 for the 2D cluster state. 



16 with a fidelity 



F 



(18) 



for N qubits, which gives an upper bound of the maxi- 
mum number of cluster qubits with a fixed fidelity 



M„ax(VV') = {l-F) 



( y 



(19) 



for V the absolute value of the average energy off- 
set of ddot pairs. In other words, approximating the 
anisotropic evolution by isotropic Ising-like evolution is 
valid provided that the total number of qubits does not 
exceed A^max, which depends on the acceptable less-than- 
unity fidelity F. 

With the evolution time t — 7r/4<;, the fidelity and the 
maximum number of cluster qubits can be written as = 
1 - (ttNT/vY and N^,^{N^) = {1 - F) [V/TiTf. For 
T = O.l^eV, = lO/xeV and T/ w 4<r = 40mcV (for large 
N), our calculations show that a 2D 162-qubit cluster 
state with a high fidelity F = 0.99 can be produced. For 



VI. CONCLUSIONS 

We have considered deterministic unitary evolution of 
a cluster state in a charge-qubit structure with charge 
qubits made of double-dot (ddot) structures. Although 
periodic evolution into charge-qubit cluster states has 
been considered before, anisotropy presents a critical yet 
overlooked challenge. At first anisotropy seems to destroy 
the opportunity to create cluster states in this way. 

We have shown how to circumvent this problem by 
applying electric field biases to the ddot structures. In 
the slow-tunneling regime, the effective single-quantized 
multi-qubit Hamiltonian can be approximated by the 
Ising-like Hamiltonian. In this case, electric field biases 
can overcome the challenge of anisotropy. If the elec- 
tric fields had to be tailored to each ddot charge qubit, 
or if the field had to be controlled over time, the strat- 
egy would be impractical. However, we have shown that 
a global field over all but the boundary qubits, and five 
choices of electric field biases on boundary qubits no mat- 
ter how large the system is. entirely eliminates the prob- 
lem of anisotropy. 

Remarkably, by changing the structure of the array of 
charge qubits, we can generate a large two-dimensional 
cylindrical cluster state by simply applying an electric 
field on the first and last rows, and a different one for 
the remaining rows, respectively. Compared to previous 
schemes, no assumption of isotropy for charge qubit cou- 
plings is made in our procedure and in fact is shown not 
to be valid for two-dimensional cluster states. We aug- 
ment our theoretical analysis of anisotropy by including 
numerical analysis for the case of GaAs double dots in a 
two-dimensional lattice. In particular we show that the 
energy offsets due to anisotropy are noneligible in this 
case. 

For these charge qubits to be useful, noise and deco- 
herence need to be considered. Also the charge qubits 
considered here periodically evolve into cluster states and 
then back to their initial states due to the periodic nature 
of the unitary evolution. Timing becomes critical in such 
dynamics, or else the interactions that produce the clus- 
ter states must be able to be switched off. Measurement- 
based quantum computing also becomes challenging for 
such periodically-occurring cluster states. These consid- 
erations are the seeds for future study. 
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